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A  nnultigrid  acceleration  technique  developed  for  solving  the  three/dimensional  Navier- 
Stokes  equations  for  subso^ic/transonic  flows  has  been  extended  to  supersonic/hypersonic  flows. 

An  explicit  multistage  Rung\rKutta  type  of  time-stepping  scheme  is  used  as  the  basic  algorithm 
in  conjunction  with  the  multigrid  scheme.  Solutions  have  been  obtained  for  a  blunt  conical 
frustum  at  Mach  6  to  demonstrate  the  applicability  of  the  multigrid  scheme  to  high-speed  flowsr 
Computations  have  also  been  performed  for  a  generic  High-Speed  Civil  Transport  configuration 
designed  to  cruise  at  Mach  3.  These  solutions  demonstrate  both  the  efficiency  and  accuracy  of  the 
present  scheme  for  computing  high-speed  viscous  flows  over  configurations  of  practical  interest.  ^ 


'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract  No. 
NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  Engineering(ICASE), 
NASA  Langley  Research  Center,  Hampton,  VA  23665. 


1 


Introduction 


During  the  last  decade  or  so,  significant  prog'^ss  has  been  made  in  the  field  of  compu¬ 
tational  fluid  dynamics  (CFD)  to  have  an  impact  on  the  design  and  analysis  of  aerodynamic 
configurations.  Solutions  of  the  Euler  (in viscid)  equations  for  essentially  complete  aircraft  con¬ 
figurations  [1-4]  and  the  solutions  of  the  Navier-Stokes  equations  for  high  Reynolds-number, 
viscous,  transonic  flows  over  aircraft  components  are  now  available  in  the  open  literature  [5-9J. 
It  is  noteworthy  thai  most  of  the  efficient  numerical  schemes  for  solving  aerodynamic  flows 
rely  on  multigrid  acceleration  technique  [2,8,9]  to  enhance  the  convergence  rate.  The  multigrid- 
based  schemes  have  the  desirable  property  that  the  number  of  iterations  required  to  achieve  a 
steady-state  solution  is  nearly  independent  of  the  mesh  size  for  a  given  class  of  problems.  Thus 
one  can  achieve  essentially  grid-converged,  steady-state  solutions,  even  for  the  numerically  de¬ 
manding  problem  of  high  Reynolds-number,  transonic,  viscous  flow  over  realistic  aerodynamic 
configurations  with  a  reasonable  amount  of  computer  resources  [9]. 

Despite  the  progress  achieved  in  solving  transonic  flows,  the  development  of  CFD  methods 
for  supersonic/hypersonic  flows  seems  to  be  lagging  behind  at  the  present  time.  With  the  current 
interest  in  high-speed  vehicles  such  as  the  High-Speed  Civil  Transport  (HSCT)  and  the  National 
Aero-space  Plane  (NASP),  it  is  imperative  that  efficient  computational  algorithms  be  developed 
for  high-speed  flow  regimes.  In  this  paper,  we  will  discuss  the  progress  made  in  this  direction 
using  a  multigrid-based  central-difference  scheme  and  present  sample  results  for  problems  of 
practical  interest. 

Governing  Equations  and  Numerical  Method 

The  basic  equations  under  consideration  here  are  the  unsteady  Navier-Stokes  equations. 
These  are  specialized  to  a  body-fitted  coordinate  system  (^,  where  //,  and  c,"  represent  the 
streamwise,  normal,  and  spanwise  coordinates,  respectively.  The  //  coordinate  lines  are  nearly 
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orthogonal  to  the  wing  surface.  Since  the  dominant  viscous  effects  for  high-Reynolds-number 
turbulent  flows  arise  from  viscous  diffusion  normal  to  the  body  surface,  a  thin-layer  assumption  is 
employed  here  by  retaining  only  the  viscous  diffusion  terms  in  the  direction.  These  equations 
can  be  written  in  the  conservation  law  form  as: 


h  + 1 


dG  dH 

-L  -  -f  - 

di)  5C 


dG„ 

dr] 


(1) 


where  the  dependent  variable  vector  U  is  given  by  the  relation 


U  =  { 


‘  P 


pu 

pv  ^ 
pro 

U-BJ 


(2) 


In  the  eq.  (1),  F,  G,  G'v,  and  H  are  the  flux  vectors,  and  J  is  the  Jacobian  of  the  transformation. 
The  complete  forms  of  these  quantities  are  readily  available  in  [9]. 

A  pseudo  time-stepping  scheme  based  on  a  Runge-Kutra  scheme  [10,11]  is  used  for  integrate 
ing  the  time-dependent  equations  to  steady  state  and  as  a  smoother  in  the  multigrid  scheme.  For 
convenience,  let  us  first  write  the  discretized  form  of  the  governing  equations  in  the  following 
operator  notation; 


^  +  Q(f)  -  DW)  =  0 


(3) 


where  Q  contains  all  the  convective  and  viscous  fluxes  and  D  represents  the  artificial  dissipative 
fluxes. 

bince  our  primary  interest  here  is  to  obtain  solutions  for  viscous  flows  via  the  Navier-Stokes 
equations,  both  the  diffusion  and  the  convective  terms  are  important  in  contrast  with  the  Euler 
equations  where  con  .'ective  terms  are  dominant.  Therefore,  it  is  preferable  to  employ  a  scheme 
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that  has  a  larger  stability  bound  along  the  real  axis  in  addition  to  good  stability  properties 
along  the  imaginary  axis.  Based  on  the  Fourier  stability  analysis  of  a  one-dimensional  model 
problem,  the  five-stage  Runge-Kutta  scheme  proposed  by  Jameson  [12],  with  3  evaluations  of 
the  dissipative  operator  at  the  first,  third,  and  fifth  stages,  appears  to  be  very  attractive  and  is 
employed  in  the  present  work.  The  convergence  to  steady  state  is  enhanced  via  the  use  of  local 
time-stepping  and  implicit  residual  smoothing  techniques  [10,11],  with  the  coefficients  of  the 
residual  smoothing  computed  in  the  manner  described  by  Vatsa  and  Wedan  [9]. 

Artificial  Dissipation  Model 


The  basic  artificial  dissipation  model  used  in  this  study  is  patterned  after  the  work  of  Jameson, 
Schmidt  and  Turkel  [10]  and  of  Jameson  and  Baker  [11]  for  2-D  and  3-D  Euler  equations 
respectively,  and  modified  later  by  Vatsa  and  Wedan  [9]  for  3-D  Navier-Stokes  equations.  In 
order  to  discuss  the  modifications  required  for  supersonic/hypersonic  flows,  let  us  first  examine 
the  dissipation  terms  in  the  i-direction: 


*+2  0'^' 


(4) 


In  the  above  expression,  1.  ^  ^  is  the  modified  eigenvalue  scaling  factor  [9]  and  the  coefficients 
and  are  related  to  the  pressure  gradient  as  follows: 

I  I  2 


(5) 


where  the  coefficients  and  are  set  equal  to  1/2  and  1/64,  respectively.  The  term  1/ 
depends  on  the  pressure  gradient  and  is  modified  to  give  a  TVD  variation  of  the  shock  switch 
[13]  in  the  following  manner: 


!/,  = 


1 7^1+ 1. j. A'  -Pi.j.k  "F  Pt—\.j.k 


( 1  -  -  Pi-i.j.aI)  +  <^{Pt+l,j,k  +  '^Pi.j.k  +  Pi-l,j,k) 


(6) 
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Note  that  by  setting  u;  =  1,  we  can  recover  the  shock  switch  that  has  been  used  in  earlier  studies 
for  computing  transonic  flows  [9-12].  For  supersonic  and  hypersonic  flows,  where  shocks  are 


much  stronger,  we  use  vj  =  1/2.  The  expressions  for  the  dissipation  terms  in  the  j  and  k 
directions  are  derived  in  a  similar  manner. 

Evaluation  of  Time-Step 


It  is  very  important  to  estimate  the  allowable  time-step  as  accurately  as  possible  in  order 
to  construct  a  robust  time-stepping  scheme.  Failure  to  do  so  generally  creates  difficulties  when 
the  scheme  is  applied  to  different  flow  problems  with  widely  varying  test  conditions  and  grid- 
densities.  An  attempt  is  made  here  to  derive  the  expressions  for  allowable  time-step  for  the 
present  scheme  from  stability  considerations.  For  convenience,  let  us  start  with  the  Navier- 
Stokes  equations  written  in  non-conservative  form: 


dU  JU  JU  M 

at  dx  dy  cjz  dx-  dy-  dz- 

,  d-^-U  dH^ 

where  U  is  the  velocity  vector  and  A,  B, - Gzx  are  the  coefficient  matrices,  their  full  form 

being  available  in  [14  ]. 

Transforming  these  equations  to  the  body-fitted  curvilinear  coordinate  system  tj  and 
and  making  the  thin-layer  assumption,  one  gets 

+  B^y  +  +  [Aiijt  +  Bijy  + 

+  [AQ  +  BQ  +  CX-]-^  (8) 

o  ■}  >  d~u 

=  [Dy'j.  -1-  Erjy  -f  Fl}Z  +  Crylh^y  +  Gy:1]y)]:  -f  GxitJ^TJ^] 

In  general  it  is  very  difficult  to  derive  an  exact  expression  for  the  time-step  At,  since  the 
coefficient  matrices  are  non-symmetric.  Abarbanel  and  Gottlieb  [14]  have  recently  delineated  a 


procedure  to  symmetrize  all  of  the.ie  coefficient  matrices  simultaneously.  Taking  advantage  of 


the  work  of  [14],  and  the  property  that  the  norm  of  a  symmetric  .natrix  equals  its  spectral  radius, 
one  can  find  an  upper  bound  on  At  in  the  foilowing  manner: 


AZ  ~  AZ^c  AZ,,  AZ(^  ''  ^ 

where,  the  first  3  terms  on  the  right  hand  side  arise  due  to  the  convective  terms  and  the  last  term 

is  due  to  diffusion  terms.  Setting  the  bounds  of  these  components  of  At  equal  to  their  respective 
spectral  radii,  we  arrive  at  the  following  expressions  for  the  convective  terms: 


=  |«^j.  +  v^y  +  +  c^J +  (f: 

^  >  A,,  =  +  vtiy  +  m//.|  +  c^7f^  +  7]}j  +  7]l 

>  A^'  =  |h(,V  +  I’Cj,  +  +  CyJCx  +  (,j  +  Ci 


(10) 

(11) 

(12) 


where  c  is  the  speed  of  sound.  The  diffusion  limit  on  the  time  step  is  obtained  in  a  similar 


manner  and  is  given  by  the  relation: 
^  Aj/// 


^Uiff 


=  ^  +  vl  +  '/fz  I  I  +  +  Iv/.-i  +  \rhnA] 


(13) 


For  viscous  flow  problems,  the  most  restrictive  time  step  is  in  the  boundary  layer  region  near 
a  solid  surface,  wnere  At,^  and  At(i,ff  are  the  most  critical  terms  in  determining  the  actual  value 
of  At.  For  adiabatic,  high  Mach  number  flow  near  a  stagnation  point,  one  can  show  that  [13] 

A,/,//  ^  A, ,.1/4  (14) 

Thus  the  diffusion  limit  on  time  step  becomes  increasingly  more  imponant  as  the  flow  Mach 
number  is  increased.  As  expected,  the  diffusion  limit  reduces  the  allowable  time  step  in  the 
near  wall  region.  However,  by  incorporating  the  diffusion  limit  in  the  time  step,  we  were  able 
to  perform  the  computations  over  a  large  range  of  Mach  number  without  changing  the  CFL 
number,  which  is  very  desirable  from  a  practical  point  of  view  since  it  helps  to  make  the  scheme 


more  robust. 


Multigrid  Acceleration  Technique 


The  convergence  acceleration  due  to  the  Ube  of  multigrid  techniques  has  been  demonstrated 
for  both  inviscid  and  viscous  flows  [8,9,12]  in  the  transonic  flow  regimes.  In  the  current 
application,  the  Full  Approximation  Storage  fFAS)  scheme  of  Brandt  [15]  is  used  in  conjunction 
with  the  multigrid  strategy  devised  by  Jameson  [12]  for  the  solution  of  the  Euler  equations. 
The  extension  of  the  scheme  of  [12]  to  the  three-dimensional  thin-layer  Navier-Stokes  equations 
for  transonic  flows  described  by  Vatsa  and  Wedan  [9]  is  used  as  a  starting  point  for  this  work. 
A  5-stage  Runge-Kutta  scheme  with  coefficients  selected  to  provide  optimum  damping  of  the 
high  frequency  errors  is  employed.  The  restriction  operator  used  to  transfer  the  solution  to  a 
coarser  grid  is  a  volume-weighted  average  of  the  eight  surrounding  cell-centered  values.  The 
forcing  function  for  a  cell  on  the  coarse  grid  is  obtained  by  simply  summing  the  residuals  of 
its  constituent  fine-grid  cells.  The  corrections  are  transferred  from  the  coarse  grid  back  to  the 
fine  grid  (or  prolonged)  by  simple  trilinear  interpolation  in  computational  space.  On  highly 
stretched  or  nonuniform  grids,  this  prolongation  operator  can  introduce  high-frequency  errors 
back  to  the  fine  grid,  causing  degradation  of  the  convergence  rate.  To  prevent  this,  the  coarse- 
grid  corrections  were  processed  through  an  implicit  residual  smoothing  operator  before  adding 
to  the  fine-grid  corrections.  Whereas  the  smoothing  of  the  coarse-grid  corrections  was  certainly 
helpful  for  transonic  flow  calculations,  it  was  found  to  be  essential  for  obtaining  converged 
solutions  for  higher  speed  flows. 

The  solutions  presented  in  this  paper  were  obtained  using  a  W-cycle,  in  which  governing 
equations  are  solved  only  in  the  restriction  step.  The  W-cycle  resulted  in  approximately  a 
25-percent  improvement  in  computational  time  compared  to  a  standard  V-cycle  for  achieving 
comparable  convergence  levels  of  the  residuals.  In  addition,  global  properties  such  as  lift  and 
drag,  develop  more  rapidly  with  the  W-cycle,  since  more  time  is  spent  on  the  coarser  grids.  It 
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was  also  found  helpful  to  run  more  cycles  on  coarser  grid  levels  for  supersonic  and  hypersonic 
flows  in  order  to  better  establish  and  precondition  the  flow  field  before  staning  computations  on 
the  finest  mesh  in  the  Full  Multigrid  (FMG)  cycle. 

The  variable-coefficient  residual  smoothings  were  applied  on  all  grid  levels  of  the  multigrid 
cycle.  On  the  finest  grid,  the  blend  of  second-  and  fourth-difference  artificial  dissipation  discussed 
previously  was  employed.  For  the  coarser  grids,  a  fixed  coefficient  second-difference  dissipation 
model  has  been  found  adequate  for  transonic  flow  computations  [9,12].  However,  the  coarse-grid 
dissipation  model  had  to  be  modified  via  a  pressure  gradient  based  TVD  switch  (see  eqs.  (5-6))  to 
improve  the  convergence  rate  of  the  present  scheme  in  supersonic  and  hypersonic  flow  regimes. 

Results  and  Discussion 

Two  test  cases  covering  supersonic  to  low  hjperscaic  speed  regimes  are  chosen  for  testing 
the  multigrid  Navier-Stokes  code  described  in  the  preceding  paragraphs.  The  accuracy  of  the 
computed  solutions  is  assessed  via  comparisons  with  available  experimental  data.  In  the  present 
investigation,  C-0  type  grids  are  employed.  The  computational  grids  are  generated  so  as  to 
cluster  points  in  the  appropriate  regions  to  resolve  sharp  gradients  present  therein.  In  addition, 
significant  grid  clustering  is  used  in  the  thin  region  adjacent  to  the  solid  surface  in  order  to 
resolve  the  thin  shear  layers  present  in  high  Reynolds-number  turbulent  flows. 

a.  Conical  Frustum  Entry  Configuration 

The  first  test  case  used  for  this  study  is  a  simple  aerodynamic  sliape  designed  for  entry  at 
Mach  5.  The  configuration  consists  of  a  modified  conical-frustum  [16]  and  is  shown  in  Fig. 
(1),  along  with  tne  downstream  and  s>mmetry  planes.  This  configuration  evolved  through  a 
conceptual  study  for  the  design  of  a  vehicle  to  accommodate  an  8-person  crew,  which  could 
sustain  a  supersonic/hypensonic  reentry  and  be  capable  of  landing  as  a  paraglider  [16].  A  wind 
tunnel  model  was  tested  at  .U^  —  (i.O  and  Reynolds  number  of  0.8x10^  based  on  the  model 
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length.  The  Navier-Stokes  calculations  were  performed  on  a  grid  consisting  of  161x65x29  mesh 
points.  The  lift-  and  drag-coetficient  data  from  the  wind-tunnel  tests  are  available  for  up  to 
an  angle  of  attack  (o)  of  12®.  The  Navier-Stokes  solutions  spanning  this  entire  angle-of-attack 
range  have  been  obtained  to  assess  the  performance  of  the  current  scheme  over  such  a  large 
range  of  test  conditions. 

The  convergence  history  in  terms  of  the  residual  error  of  the  continuity  equation  and  the  lift 
coefficient,  C/  for  the  q  =  6®  case  are  shown  in  Fig.  (2)  as  a  function  of  work-units,  where 
a  work  unit  represents  the  computational  effort  required  for  one  fine-mesh  iteration.  A  total  of 
400  iterations  (620  work-units)  were  performed  on  the  fine  grid  which  resulted  in  approximately 
seven  orders  of  reduction  in  the  residual.  The  lift  and  drag  coefficients  for  this  case  converged 
to  within  0.1%  of  their  final  values  in  less  than  50  fine-grid  iterations.  The  convergence  histories 
for  this  series  of  test  cases  are  similar  to  the  o  =  6®  case  shown  in  Fig.  (2)  except  for  the  case 
of  Q  =  12®,  for  which  the  residual  started  oscillating  after  dropping  approximately  five  orders, 
possibly  due  to  slight  unsteadiness. 

The  computed  lift  and  drag  coefficients  for  a  =  0-12®  shown  in  Fig.(3),  are  found  to  be  in 
excellent  agreement  with  the  experimental  data  over  the  complete  range  of  a.  Thus,  not  only  does 
the  present  scheme  have  good  convergence  properties  for  these  test  conditions,  but  in  addition  it 
produces  accurate  solutions  that  are  in  good  quantitative  agreement  with  the  experimental  data. 

A  better  understanding  of  the  overall  flow  field  is  obtained  by  examining  Fig.  (4),  where 
the  pressure  contours  for  the  a  =  6®  case  are  shown  on  the  symmetry  and  downstream  planes, 
in  addition  to  the  body  surface.  One  can  clearly  observe  the  nearly  conical  growth  of  the  shock 
surface  in  the  streamwise  direction.  In  the  crossflow  direction,  the  distance  between  the  shock 
and  the  body  surface  goes  through  a  minimum  near  the  tip  and  then  increases  towards  both  the 
windward  and  the  leeward  planes  of  symmetry. 
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b.  High-Speed  Civi'  Transport 

The  next  test  case  considered  here  is  that  of  the  flow  over  a  generic  High-Speed  Civil 
Transport  (HSCT)  configuration  designed  to  cruise  at  a  Mach  number,  =  3.0.  The 
conceptual  development  and  geometric  details  of  this  highly  blended  wing/body  configuration 
are  available  in  [17].  A  wind-tunnel  model  representative  of  this  vehicle  was  tested  over  a 
large  angle-of-attack  range  for  several  Mach  numbers,  and  the  experimental  data  from  this  study 
have  been  documented  by  Covell  et  al.  [18].  In  the  present  study,  we  will  concentrate  on 
the  =  3.0  case,  which  was  the  design  cruise  Mach  number  for  this  configuration.  The 
corresponding  Reynolds  number  based  on  model  length  was  6.3x10^. 

A  grid  consisting  of  145x65.x73  nodes  was  generated  for  the  Navier-Slok.es  calculations.  A 
partial  view  of  the  mesh  on  the  symmetry  plane  and  several  stream  wise  cuts  is  displayed  in 
Fig.  (5).  The  outer  boundaries  of  the  g.id  were  placed  so  as  to  contain  the  shock  emanating 
from  the  leading-edge  within  the  computational  domain.  Grid  clustering  in  the  tip,  leading  and 
trailing-edges  are  used  to  accurately  resolve  the  flow  in  high  gradient  regions. 

The  convergence  histories  for  the  lift  and  residual  of  the  continuity  equation  for  the  a  =  -5“ 
are  shown  in  Fig.  (6),  which  shows  a  reduction  of  three  orders  in  the  residual  after  300  fine 
grid  iterations  (470  work-units).  The  lift  for  this  case  converged  after  about  50  iterations.  The 
convergence  histories  for  the  entire  angle-of-attack  range  considered  here  (a  =  0  —  b")  are  very 
similar  to  this  case.  The  integrated  lift  and  drag  »,oefficients  are  con.pared  with  the  experimental 
data  of  [18]  over  the  complete  range  of  angle  of  attack  in  Fig.  (7).  It  is  clear  from  this 
figure  that  the  computed  results  are  in  very  good  agreement  with  the  experimental  data.  This  is 
very  encouraging  especially  for  drag  prediction,  since  viscous  drag,  which  is  difficult  to  predict 
accurately,  constitutes  a  significant  part  of  the  total  drag  for  this  configuration  at  Mach  3. 

Next  we  examine  the  detailed  surface  pressure  distributions  for  this  vehicle.  For  this  purpose. 
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we  concentrate  again  on  the  a  =  o"  case  for  which  experimental  pressure  data  are  available  at 
selected  streamwise  and  spanwise  stations.  Since  the  computational  grid  does  not  follow  the 
cuts  along  which  experimental  data  was  acquired,  the  computed  pressure  distributions  had  to 
be  interpolated  for  a  meaningful  comparison.  The  interpolated  solutions  at  two  x-locations  are 
compared  with  the  experimental  data  in  Fig.  (8).  For  the  station  x=129,  experimental  pressures 
are  available  only  on  the  upper  surface,  whereas  both  upper-  and  low'er-surface  pressures  are 
available  at  the  x=171  station.  The  agreement  with  the  measured  pressure  data  at  both  of  these 
stations  is  quite  good.  It  is  noted  that  the  computed  solutions  predict  the  correct  variation 
in  pressure  even  in  the  vortex-dominated  flow  near  the  wing-tip  region.  Similar  comparisons 
have  also  been  obtained  at  o  =  1“  and  «  =  T';  however  these  results  are  not  shown  here  for 
conciseness. 

Whereas  the  global  force  coefficients  and  surface  pressure  comparisons  are  helpful  in  quan¬ 
titative  validation  of  a  prediction  method,  these  flow  properties  are  inadequate  for  understanding 
the  true  three-dimensional  nature  of  the  flow  problem,  such  as  the  development  of  vortical-flow 
regions  off  the  wing-tips.  Following  the  lead  of  119],  an  attempt  is  made  here  to  visualize  the 
vortical  flow  by  plotting  the  density  contours  at  fi.xed  x-locations.  This  is  done  in  Fig.  (9)  for 
x=129  and  x=I71,  the  same  stations  for  which  the  surface  pressure  distributions  were  exam¬ 
ined.  The  experimental  laser  sheet  photographs,  that  are  used  routinely  for  visualizing  shock 
and  vortex  formations  are  also  shown  in  Fig.  (9)  at  the  corresponding  stations.  The  core  of 
the  vortex  forming  between  the  wing-tip  and  the  fuselage  and  its  feeding  sheet  are  the  most 
dominant  flow  sumeture  visible  in  these  figures.  After  a  careful  examination  of  these  results, 
the  following  observations  can  be  made.  First,  the  core  of  the  vortex  is  predicted  to  grow  in 
size  and  lift  off  farther  away  from  the  surface  as  one  moves  downstream.  Secondly,  the  overall 
shape  and  size  of  the  vortex  predicted  by  the  Nav  icr-Stokes  calculations  are  in  reasonably  good 
agreement  with  the  la.ser-sheet  data. 
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Concluding  Remarks 


A  multigrid  acceleration  technique  developed  originally  for  transonic  flows  has  been  extended 
to  solve  the  three-dimensional  Navier-Stokes  equations  for  supersonic  and  hypersonic  viscous 
flows.  The  convergence  rate  of  the  modified  mu'tigrid  code  for  obtaining  steady-state  solutions 
of  high  Reynolds-number  viscous  flows  in  the  supeisoniL/hypersonic  flow  regimes  (  up  to  Mach 
6),  is  comparable  to  the  convergence  rate  in  transonic  flows.  Convergence  histories  and  detailed 
comparisons  with  the  experimental  data  are  presented  for  two  problems  of  practical  interest. 
Based  on  these  solutions,  it  is  concluded  that  the  resulting  code  is  capable  of  predicting  high¬ 
speed  viscous  flow  problems  in  an  efficient  and  reliable  manner. 
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Fig.  2  :  Convergence  history  for  conical  frustum  on  161x65x29  grid, 

Moo  =  6.0,0  =  =  0.8x10*^ 
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Fig.  5  :  Partial  view  of  grid  for  IISCT 
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Fig.  6  :  Convergence  h-story  for  MSCT  on  145x65x73  grid, 
yUoo  =  3.0,  a  =  5",/?c/  =  6.3x1 0'^ 
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Pig.  9  :  Flow-field  visualization  for  HSCT, 
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